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Abstract: 

We highlight an interesting mapping between the moving breather solutions of the generalized Nonlinear 
Schrodinger (NLS) equations and the static solutions of neutral scalar field theories. Using this connection, 
we then obtain several new moving breather solutions of the cubic-quintic NLS equation both with and 
without uniform phase in space. The stability of some stationary solutions is investigated numerically and 
the results confirmed via dynamical evolution. 
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1 Introduction 



The nonlinear Schrodinger equation is one of the most celebrated integrable nonlinear equation which has 
found application in several areas of physics including the self-focusing of intense Laser beams, Langmuir 
waves in Plasma, the collapse of the Bose-Einstein condensates (BECs) with attractive interactions etc 
pp. While most of the applications so far are those of NLS with cubic nonlinearity, in recent years the 
cubic-quintic NLS (CQNLS) has also found several applications. Some of these are in nonlinear optics 
and BEC. For example, in nonlinear optics, it describes the propagation of pulses in double-doped optical 
fibres [2] and in Bragg-gratings [3]; in BEC it models the condensate with two- and three-body interactions 
[11[3]. It is thus of considerable interest to obtain exact solutions of the CQNLS and other generalized NLS 
equations. 

In particular, a considerable amount of attention has been given in the recent decades to localized 
solutions of various NLS, although these solutions exist as the hyperbolic limit of appropriate general elliptic 
function solutions. This direction has been explored in the context of an external periodic potential [6] as 
well as weakly interacting solitary waves in generalized NLS equations [?]• Herein we explore exhaustively 
the analytic elliptic function solutions in the absence of an external potential. 

In this context it is worth noting the connection between the moving breather solutions of the generalized 
NLS equations and the (static) solutions of neutral scalar field theories, known as Galilean invariance. 
Consider the generalized NLS equation 

iut + Uzz + f{\u\'^)u = 0, (1) 
where /(.) is a real valued algebraic function with /(O) = 0. On using the ansatz 

u{z, t) = (j){z -vt- zo) exp{{i/2)v{z - kt - 5)} , (2) 
where k = (w^ + Aa)/2v, it is easily shown that u satisfies Eq. ([1]) provided (p{x) satisfies the equation 

•/"xx + + a0 = 0, x = z-vt-zo, (3) 

which is the static field equation for a scalar field theory with nonlinear interaction term f[(/)'^)(/). This 
is a travelling wave solution, and can also be posed in x coordinates by the transformation dz — )• dx and 

2 



dt — )• —vdx + dt, in which case, considering zq = 0,5 = without any loss of generahty, Eq. ([2]) can be 
considered separable in this frame as 

u{x, t) = (/>(x) exp {{i/2)vx} exp {(i/2)(u^ — 4a)f} = exp {— mt} . (4) 

Now this Galilean transformed solution is a stationary solution in the traveling frame. Actually, we 
can slightly generalize the above ansatz, i.e. consider instead 

u{z, t) = (j){z — vt — zq) exp {{i/2){vz — vkt — 5 + r]{x))} , x = z — vt — zq . (5) 

On equating real and imaginary parts, it is easily shown that, given g, if 

dx' 

and 

(t>xx + f{(t>^)(t> + a4>= ^ , x = z-vt- ZQ, (7) 

with 4a = v{2k — v) as above, then ([7]) solves ([T]). (Note in the traveling frame of x there would also be an 
extra nonlinear term due to the first derivative, for v ^ 0.) 

It thus follows from here that knowing the various static solutions of the scalar field theory ([3]) or 
([7|), one can immediately write down the moving breather solutions of the corresponding NLS model by 
using Eqs. ([T]) and ([2]). Note that the breather solution of the corresponding NLS model will be uniquely 
characterized by its five parameters: v,a,g,6 and zq. In particular it follows from here that in order to 
obtain the breather solutions of the CQNLS equation, which is characterized by 

f{\u\^) = b\u\' + c\u\\ (8) 

one needs to obtain the static solutions of the (fp'-(f)'^-(f)^ field theory characterized by 

(t>xx + a(l) + h<f + = , (9) 

in case 5 = while if 5 7^ 0, then one needs to obtain the solutions of the equation 

(i>xx + a(l) + h4>'^ + c4f' - -^ = 0. (10) 



Note that in the special case of c = 0, one gets the well known mapping between the celebrated (cubic) 
NLS and (p'^-cp^ field theory. If instead one wants to obtain the breather solutions of the quadratic-cubic 
NLS characterized by 

f{\u\^) = b\u\+c\u\\ (11) 
then one needs to obtain the static solutions of the (fp'-cp'^-cf)'^ field theory characterized by 

(t>xx + a(t) + bcp"^ + c(p^ = 0. (12) 

We have thus transformed the problem of finding the breather solutions of the generalized NLS models 
to that of finding static solutions of neutral scalar field theories with power law potentials. There is 
however, one important difference. In view of the fact that the potential must be bounded from below, in 
neutral scalar field theory models, normally one takes the coefficient of the leading term in the potential 
to be positive, i.e. in Eqs. ([9|) and (jl2p . one takes the coefficient c < 0. However, in the context of the 
generalized NLS equations, this coefficient need not necessarily be negative. In fact usually it is taken 
to be positive, yielding bright solitons while if it is negative, then one obtains dark solitons. It is easily 
seen that taking the coefficient c > in Eqs. Q and p2|) is equivalent to considering the solutions of the 
nonlinear oscillator problem (rather than the field theory problem). 

The purpose of this paper is to give exhaustive solutions to the CQNLS as well as quadratic-cubic 
NLS equation as given by Eq. (|12|) . For that purpose we make use of the known static solutions of the 
corresponding cp'^-cp'^-cj)^ as well as (fp'-cp'^-cf)'^ field theory as well as oscillator models and also obtain few 
new solutions, which to the best of our knowledge, had not been explicitly written down before in the 
literature. Further, we also obtain static solutions of the field theory c/i^-c/i^-^^-^"^. In this context, it is 
worth pointing out that recently we have obtained several static solutions of the coupled (jp'-cp'^-cp^ [8j, <fP'- 
(p^-cj)^ [9] as well as (fp'-<j)^ |10j models from where we can immediately obtain the solutions in the decoupled 
case, for c < 0. Extending these ideas, we also obtain solutions of these models even in case c > 0. 

The paper is organized as follows. In Sec. [2] we obtain static solutions of the constant (or linear in the 
case V 7^ 0) phase CQNLS problem (i.e. when (/>~^ term is absent). In section [2^2] we investigate the linear 
stability of some of these solutions numerically. In particular, we focus on the hyperbolic limit in which 
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our results here are consistent with the well-known Vakhitov-Kolokolov criterion [11] as well as earlier 
results for a more general non-linearity [121 US]- Dynamical evolution confirms the stability results and a 
connection is observed between the nature of the nonlinearity and the behavior of the unstable evolution. 
In particular, in the case c > 0, when 6 < the solutions tend to blow-up with a self-similar behavior, 
while unstable solutions for which 6 > do not (perturbations are additive). Then, in Sec. Owe provide 
solutions for the (j)"^ — (p'^ — (p^ — model as given by Eq. Q in case each of a, b, c can be either positive 
or negative. In Appendix [X] we provide static solutions for the cjp' — (p"^ — cp'^ model as given by Eq. p2|) in 
case each of a, b, c can be either positive or negative. These will be relevant in the context of the breather 
solutions of quadratic-cubic NLS. For completeness, in Appendix [B] we also provide static solutions for the 
0^ — (p'^ model as given by Eq. (jlip in case each of a, c can be either positive or negative. These will be 
relevant in the context of the breather solutions of NLS. Finally, in Sec. |l]we summarize the results and 
indicate possible relevance of these results. 



2.1 Theoretical development 

Let us consider solutions of field Eq. ([9]). As explained above, once these solutions are obtained, then the 
solution of the CQNLS equation are immediately obtained from here by using Eqs. ([T]), ([2]) and ([8]). We 
list below eight distinct periodic solutions to the field Eq. ([9]), i.e. the trivial phase case of g = 0. In each 
case, we also mention the values of the parameters a, b, c, in particular, if they are positive or negative. 

2.1.1 Dark soliton families 

Solution I 

It is easily shown [8] that 



2 Solutions of 02_04_06 hence CQNLS Model 



(p = A\/ 1 -|- sn(i?x -t- xo, m) 



(13) 



is an exact solution to the field Eq. Q provided 



(5m - 1)5^ = -4a, {hm-\)J^ 



8ma 



3{5m — 1)6^ = 64mac . 



(14) 



b 
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Thus this solution is vahd provided 6 > 0, c < while a > (< or =) depending on if m < (> or =) 1/5. 
Here sn(a;,m) and cn(x,m), dn(3;,m) denote Jacobi elliptic functions with modulus m [14:\ I15|. 
In the limit m = 1 the periodic solution (fTBj) goes over to the dark soliton solution 

(p = Ay/1 + tanh{Bx + xq) , (15) 

provided 

B'^ = -a, A'^ = -{2a/b), b"^ = {16/3)ac . (16) 

Thus the dark soliton solution exists to field Eq. ([9]) provided a < 0, 6 > 0, c < 0. 
Solution II 

It is easily shown [8j that 



Asn{Bx + xq, m) 



(17) 



Y^l — Dsn^ {Bx + xo,m) 
is an exact solution to the field Eq. ^ provided 

[3D - (1 + m)]B^ = -a , = 2[2D{1 + m)-m- ?>D'^]B'^ , 

36^ _ [2D(1 + m)-m- ZD'^f 
4^ ~ D{1 - D){m - D)\3D - (1 + m)] ' 

There are different constraints depending on the value of D. For example if D < 0, then this solution is 
valid provided a > 0, 6 < 0, c > 0. On the other hand, if D > then while the solution is only valid if 
c < 0, the signs of a, 6 depend on the value of D. For example, while a < (> or =)0 depending on if 



(18) 



D > (< or =) (1 + m)/3, 6 < (> or =) depending on if D < (> or =) Vi m +m^^ Summarizing 

„ 1 + m — y/l — m + m?' 
0<D< , a>0,6<0,c<0. 



1 + m — \/l — m + m? 1 + m 

<D< — - — , a>0,6>0,c<0, 

o o 

2 -|- 

<D<m, a<0,6>0,c<0, 

3 

D<0, a>0,6<0,c>0. (19) 



In the limit m = 1 the periodic solution ()17p goes over to the dark soliton solution 

A tanh(Bx + xq) 
1 - Dtanh^{Bx + xq) 



(20) 



provided 



{3D-2)B^ = -a, bA^ = 2{3D - - D)B\ — = . (21) 



4ac D{3D - 2) 

Thus the dark sohton solution exists to field Eq. Q provided the following constraints are satisfied 

0<L'<^, a>0,b<0,c<0, 
1 2 

-<D<-, a>0, 6>0, c<0, 

o o 

2 

-<i:'<l, a<0,6>0,c<0, 
o 



i:»<0, a>0,6<0,c>0. (22) 



Solution III 

It is easily shown ^ that 



(23) 



(24) 



-y/l — Dsn^{Bx + xo, m) 
is an exact solution to the field Eq. ([9]) provided 

[3m - (1 + m)D]B'^ = -aD , DbA^ = 2[D^ - 2D{1 + m) + 3m]B^ , 

36^ _ [D"^ - 2D{1 + m) + 3m]2 
4^ ~ (1 - D){m - D)[3m - (1 + m)D] ' 

There are different constraints depending on the value of D. For example if I? < 0, then this solution is 

valid provided a > 0, 6 < 0, c > 0. On the other hand, if D > then while the solution is only valid if 

c < 0, the signs of a, b depend on the value of D. For example, while a < in case m > 1/2, for m < 1/2, 

a < (> or =) depending on if D < (> or =) ^j^p^- On the other hand, 6 > (< or =) depending on if 

D < {> OT =)1 + m — \/l — m + m?. 

In the limit m = 1 the periodic solution ([23]) goes over to the dark soliton solution 

A , ^ 

4> = I , (25) 

Jl- Dianh^{Bx + XQ) 

provided 

{3-2D)B' = -aD, DbA^ = 2{3 - D){1 - D)B\ _ = . (26) 

There are different constraints depending on the value of D. For example if D < 0, then this solution 
is valid provided a > 0, 6 < 0, c > 0. On the other hand, if D > then the solution is only valid if 
a < 0, 6 > 0, c < 0. 



It is worth pointing out that for D < 0, solution (j23p is not an independent solution but rather can be 
easily derived from the solution (|37|) below by using well known transformation formulas for Jacobi elliptic 
functions and making use of the translational invariance of the solutions. However, for D > 0, (j23|) is an 
independent solution. 

2.1.2 Bright soliton families 

Solution I 

It is easily shown that 

(j) = A^J\ + cn(Bx + xo, m) , (27) 
is an exact solution to the field Eq. ([9]) provided 

(4m + \)B^ = 4a , (4m + l)A^ = — , 3(4m + 1)6^ = 64mac . (28) 

Thus this solution is valid provided a > 0, 6 < 0, c > 0. 

In the limit m = 1 the periodic solution (p7|) goes over to the bright soliton solution 



(/> = Ay/1 + sech{Bx + xq) , (29) 

provided 

= 4a , = -(8a/6) , 156^ = 64ac . (30) 

Thus the bright soliton solution exists to field Eq. ([9]) provided a > 0, 6 < 0, c > 0. 
Solution II. 1 
It is easily shown that 

(f) = AyJ\ + ^u{Bx + Xo, m) , (31) 
is an exact solution to the field Eq. ([9]) provided 

(4 + m)B2 = 4a, (4 + m)^2 = -y, 3(4 + m)^^ = 64ac . (32) 

Thus this solution is valid provided a > 0, 6 < 0, c > 0. 

In the limit m = 1 the periodic solution (j3ip goes over to the bright soliton solution (j29p satisfying the 
constraints (l30l). 



Solution II.2 

It is easily shown that 

4> = Ay/dn{Bx + xo,m) + k' , k' = \/l - m , (33) 
is an exact solution to the field Eq. Q provided 

(4 - 5m)B'^ = 4a , bA^ = -2B'^k' , cA^ = (3/4)5^ . (34) 
Thus this solution is valid provided 6 < 0, c > while a > (< or =) depending on if m < (> or =) (4/5). 

In the limit m = 1 the periodic solution (|33p goes over to the bright soliton solution 

(/) = A\J sech(i?x + xo) , (35) 

provided 

= -4a , 6 = 0, cA^ = (3/4)5^ . (36) 

Thus the bright soliton solution exists to field Eq. Q provided a < 0, 6 = 0, c > 0. 
Solution II.3 
It is easily shown that 

^ = A\J dn(i?x + xo, m) — /c' , A:' = \/l — m , (37) 
is an exact solution to the field Eq. Q provided 

(4 - 5m)52 = 4a , hJ^ = 2B^k' , cA^ = (3/4)^2 . (38) 
Thus this solution is valid provided 6 > 0, c > while a > (< or =) depending on if m < (> or =) (4/5). 

In the limit m = 1 the periodic solution (j37p goes over to the bright soliton solution (j35p satisfying the 
constraints (p6]) . 
Solution III 
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It is easily shown [8] that 

Acn{Bx + xo,m) 



(39) 



(40) 



Y'T^^"Dsn2^^x~+~xo7'^ ' 
is an exact solution to the field Eq. ([9]) provided 

[2m - 1 - D{2 - m)]B'^ = -{1- D)a, (1 - D)bA^ = 2[(1 + 2D)m - D{D + 2)]B^ , 
36^ _ [(1 + 2Z))m - D(Z) + 2)]2 
4oc ~ D{m -D)[2m-1- D{2 - m)] ' 

There are different constraints depending on the value of D. For example if D < 0, then this solution is 

valid provided c > 0. On the other hand, the signs of a, h depend on the value of < 0. For example, 

while for m > 1/2, a < 0, for m < 1/2, a > (< or =) depending on if \D\ < (> or =) 2^^™ ■ On the other 

hand 6 > (< or =) provided \D\ < (> or =) 1 — m + \/l — m + . 

On the other hand, if < D < m, then this solution is valid if c < 0, while the signs of a, b depend on 

the value of D > 0. In particular, while a > in case m < 1/2, for m > 1/2, a < in case < D < ^^"^^ 

while a > provided < D < m. On the other hand, 6>OifO<L'<— (l + m) + Vl — m + , 

6 < in case — (1 + m) + \/l — m + < D < m. 

In the limit m = 1 the periodic solution (I39p goes over to the bright soliton solution 

^ Asech.{Bx + xq) ^^^^ 



'1 - Dtanh-^iBx + xo) 
provided 

B^ = -a, bA' = 2{l + D)B\ ^ = (1±^ . (42) 

Aac D 

Thus the bright soliton solution exists to field Eq. ([9|) provided the following constraints are satisfied 

0< D <1, a <0,b> 0,c <0, 
-1<D<0, a<0,6>0,c>0, 

D<-1, a<0,6<0,c>0. (43) 

Solution IV 

It is easily shown |8j that 

^_ Adn{Bx + xo,m) ^^^^ 



\/l — Dsn^{Bx + Xq, m 
10 



is an exact solution to the field Eq. Q provided 

[2m - - D{2m - 1)]B'^ = -{m - D)a , 

(m - D)bA'^ = 2[(1 - 2D)m - D{D - 2)]B'^ , 
362 ^ [(1 - 2D)m - D{D - 2)]^ 

4ac D{l-D)[2m-m? -D{2m-l)]' ^ ' 

There are different constraints depending on the value of D. For example if < D < m, then this solution 
is valid provided a < 0, & > 0, c < 0. 

On the other hand, if < 0, then this solution is valid if c > 0, while the signs of a, h depend on the 
value of < 0. In particular, while a < in case m > 1/2, for m < 1/2, a < in case < \D\ < '^l^2rn 
while a > provided ™^^2m'* < D < ni. On the other hand, 6>OifO<-D<— (1 — m) + \/l — m + 
while 6 < in case —(1 — m) + \/l — m + < D < m. 

In the limit m = 1 the periodic solution ()44p goes over to the bright soliton solution (I4ip satisfying the 
constraints as given by Eqs. and 

It is worth pointing out that for D > 0, solution ()44p is not an independent solution but rather can 
be easily derived from the solution ()39p by using well known transformation formulas for Jacobi elliptic 
functions and making use of the translational invariance of the solutions. However, for D < 0, ()44p is an 
independent solution. 

2.2 Numerical Solutions and Stability Analysis 

We will now look at some of the particular solutions presented above and investigate their stability prop- 
erties both through a systematic linear stability analysis and also dynamical evolution. In particular, we 
consider exhaustively the various regimes of hyperbolic dark and bright soliton solutions with trivial phase 
from section [21 as well as some elliptic examples from the same section. In considering the structural 
linearized stability we will consider the static reference frame (i.e. v=0) and focus on the stationary wave 
solutions. Such stationary solutions have Galilean invariance, yielding the ansatz for the traveling wave 
given in Eq. ([3]). Although, assuming a traveling reference frame and Galilean transformations of the 
perturbations as well, the equations for dynamical stability of static {v = 0) and traveling wave {v / 0) 
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Figure 1: (Color Online) Left: Dark soliton solution I with m = l,a = —2, 6 = 1 (left), and the linearization 
spectrum as computed by Eq. (|46p . Right: Dynamical evolution of n = + r, where r = pmax|(/>(x, 0)pC/, 
p = 0.05, [/ is a uniformly distributed random variable in (0, 1), and (j) is solution I with m = 1, a = —2, b = 
1, as shown to the left. The stability is confirmed as the solution sustains only minor oscillations from its 
original form after the 5% perturbation. 

solutions are equivalent. We note that above the Landau critical velocity, the solution will be so-called 
Landau-unstable [16j, although this instability is not expected to have any effect in the absence of any ex- 
ternal potential, and is usually investigated in the presence of an external impurity [17J or lattice [18]. Given 
a solution </> of ([3]), we consider the linearization of Eq. ([1]) around this stationary solution, i.e. assume 
the following ansatz, for e ^ 1, 

u{z,t) = [4>{z) + eu{z,t)]. (46) 
Assuming u = w{x)e'^^^ is separable then we have the following linearization system 





V = \v 



for any arbitrary eigenfunction v (which includes the possibility of a more general spatial form of the 
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perturbations). The blocks are given by 

A = a + 5^x + 26|(/>|2 + 3c|(/)|^ 

This eigenvalue problem, commonly known as the Bogoliubov system, has been solved with both the 
matlab functions eig in which the full matrix is diagonalized and eigs, which implements the standard 
ARPACK Arnoldi iterative algorithm to solve for the smallest/largest n eigenvalues. In particular, the 
latter, quicker and more efficient method for evaluating a particular subset of eigenvalues, is used for 
exhaustive continuations and benchmarking over much finer spatial grids, while agreement is found in all 
cases in which both methods are used. 

2.2.1 Dark soliton families 

For the first family of solutions given by Eqs. (jlSllSp . the hyperbolic tangent [m = 1, Eq. (jlSp ] (also 
known as "kink", or "dark soliton") solution was found to be stable, just like in the case of the cubic NLS. 
An example solution is presented in the left panels of Fig. [1] for a = —2 and 6=1. Continuations of 
these solutions for (a, b) S [—2, —.1] x [.01, 1] were performed and there was no deviation from this stability 
result. It seems reasonable to conclude that this is universal, although we did not do an exhaustive search 
of parameter space. The stability results were then confirmed with dynamical evolution using a standard 
fourth order Runge-Kutta method. The right panel of Figure ([T]) displays the result of evolving u = (j) + r, 
where r = pmax|</>(x, 0)p[/, p = 0.05 and C/ is a uniformly distributed random variable in (—1, 1). Notice 
even with this large perturbation, the robust structure persists for up to at least t = 1000, experiencing 
only minor oscillations around the stationary solution. 

It is most likely a result of the fact that the solutions for m ^ 1 are "very nearly non-differentiable" 
that, even with appropriately matched periodic boundary conditions, the norm difference of the exact 
solution with the one on the numerical grid is orders of magnitude larger than that of the hyperbolic kin. 
This prevents such a systematic numerical analysis of the stability, although using the exact solution as 
an initial guess for a fixed point solver on the numerical grid, the solution does converge to a smoother 
version which is wildly unstable, and this suggests these solutions are unstable. 
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Figure 2: (Color Online) The left two panels show the dark soliton solution II (left) with unstable param- 
eters {a,b, D,m) = (0.75,1,0.3,0.7) and the linearization spectrum is on the right as given by Eqs. ()46p . 
The right panels show a solution from the same family but for a set of parameters for which the solution 
is stable, {a,b,D,m) = (0.1,0.5,0.3,0.5). 
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Figure 3: (Color Online) The left panel is the dark soliton solution II with unstable parameters 
{a,b, D,m) = (0.75,1,0.3,0.7) from the left of Fig. [21 while the right is for the unstable parameters 
{a,b,D,m) = (1.5,1.5,0.5,0.99). 
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Next, we consider the smoother family of dark-sohton solutions II given in Eqs. (|17l20p . For the 
hyperbolic solutions (m=l), Eq. (|20|) . a representative value of D was chosen for all the intervals defined 
in Eq. (j22|) . with the exception of the interval 2/3 < D < 1, where the signs of b and c are the same 
as solution I above. For each of these values of D G {—0.2,0.1,0.35}, a two parameter continuation was 
done for the corresponding values of (|a|,|6|) G [-1,1] x [-01,1]. All solutions were found to be stable. 
Some periodic solutions (m ^ 1) from this smoother family were also found. For certain parameters, 
the solutions are actually stable, while for others, they are very unstable. The unstable solution for 
{a,b, D,m) = (0.75,1,0.3,0.7), and the stable one for {a,b, D,m) = (0.1,0.5,0.3,0.5) as given in Eq. 
(jl7p . and their linearization spectra are shown in Fig. ([2]). We note here that the Floquet theorem 
has been invoked in conjunction with periodic boundary conditions, since the stability matrix of the 
periodic solution is periodic, in order to scan the infinite energy spectrum with a subdivision of 200 
separate spatial frequencies out of the infinitely many. Notice the continuum of eigenvalues with negative 
imaginary part which form loops symmetric with respect to both the real and imaginary axes (the symmetry 
is a result of the fact that the matrix is Hamiltonian) . The dynamics of the unstable solution with 
{a,b, D,m) = (0.75, 1,0.3,0.7) from the left panels of Fig. [2] are shown on the left panel of Fig. [3l The 
linear instability has almost no discernible effect on the evolution even with a very large 25% perturbation 
of the initial amplitude. Somehow the stability properties of the hyperbolic limit are "inherited" by its 
family, regardless of stability of the linearized system. This presumably has to do with this solution 
residing at a small steep local maximum of the energy within a large basin, in which a small perturbation 
is enough to leave the linear regime, but not to significantly alter the structure of the solution (the unstable 
eigenfunctions are presumably of a very similar form as the solution) . We believe this "non-linear stability" 
is an interesting phenomenon, but is outside the scope of this paper and will be investigated further 
elsewhere. The right panel shows the evolution of another unstable solution from this family for parameter 
values (a, b, D, m) = (1.5, 1.5, 0.5, 0.99), where the instability has a greater effect, but still rather negligible 
and at long times considering the large perturbation. Many other unstable solutions from this family were 
evolved with large perturbations and this was the most unstable among them. 
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Figure 4: (Color Online) Left panels are the bright soliton solution I with m = 1, a = 1, 6 = —1 (left), and 
its linearization spectrum. Right panels depict the power (left, given on a log scale) and stability (right), 
as indicated by the maximum imaginary part of the spectrum, of a two-parameter family of bright-soliton 
solutions I (of Eq. (29) with m = 1 and (a, 6) G [0.4, 1] x [-0.01, -1]. 

2.2.2 Bright soliton families 

Next we turn to the bright soliton solutions. Now the hyperbolic [m = 1, Eq. [29] solution (namely "pulse", 
"bright soliton", or just "soliton") is actually unstable (see, e.g. [E]). This is markedly different from 
the cubic NLS case, where the analogous solution is stable. This can be connected to the fact that it 
is well-known that the norm (or, energy) of the solution to Eq. ([1]) for f{x) = 'J?" is bounded for 
o < 2/d, where d is the dimensionality of the problem [13]. The value a is known as the critical exponent, 
and in our case of d = 1, we have o" = 4. This is the smallest power nonlinearity for which blowup can 
occur, corresponding to an exact balance between kinetic and potential energies under the constraint of 
conserved mass. The solution and its stability are presented in the left panels of Fig. (H]). A two-parameter 
continuation in the parameters a and b (with c determined by them) was performed in order to confirm 
that the family of solutions is always unstable in the region (a, 6) G [0.4, 1] x [—1, —0.01], as presented in 
the right panels of Fig. ([!]). The stability here can also be understood in terms of the Vakhitov-Kolokolov 
criterion [llj, since dN/da > 0. The small region of stable solutions with very large amplitudes for very 
small b are actually constant solutions, since for increasing b the amplitude of the soliton shrinks, while 
the "pedestal" (plane wave) it is sitting on grows. The dynamical evolution given in Fig.dJ]) confirms the 
prediction and the solution collapses around t = 15. The collapse of the waveform is presented on a log 
scale so it can be appreciated. 
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Figure 5: (Color Online) Dynamical evolution of u = (p + r, where r = pmax|</>(x, 0)p{7, p = 0.05, [/ is a 
uniformly distributed random variable in (0,1) and (p is solution I with m = l,a = 1,6 = —1, as shown 
if Fig. The stability is confirmed and the solution collapses with exponentially self-similar nature. A 
log scale is used for the intensity to illustrate this. 
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Figure 6: (Color Online) Top: bright soliton solution II. 2 with (a, b, m) = (—0.15, —1, 0.9), (c ~ 1.56). The 
linearization spectrum is given in the middle as defined in Eqs. ()46p . and the dynamical evolution is to the 
right. Bottom: Bright soliton solution III with {a,b,m,D) = (-0.1,1,0.4,-0.3), {c^ 0.98). 
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It is interesting to note that the instabihty to collapse appears to be correlated to the sign of the 
cubic term for the bright solitons (again consistent with the findings of [12j). That is to say, those bright 
solitons we found for which 6 > were stable, while for b < they were unstable (c > of course). 
In particular, we now turn to solutions of Eq. (|39|). Two parameter continuations were performed for 
D £ {-1.5,-0.3,-0.1,0.1,0.6} and (|a|,|6|) G [0.1,1] x [0.01,1], and the last four were invariably stable, 
for all of which 6 > (for the former two c > and for the latter two c < 0). For the first value of D 
the corresponding b is negative (and of course c > 0) and, as in solution I, the solution is unstable in the 
entire parameter region. It should be noted that the solution in Eq. (|35p where 6 = and c > is linearly 
stable as the imaginary pair of eigenvalues for 6 < pass through the origin and become real for 6 < 0. 
This extra pair of eigenvalues is in addition to one pair associated to phase invariance, 6 in Eq. ([2]), and 
one pair associated to translational invariance, zq in Eq. ([2]), and is due to the conformal invariance of the 
soliton solution of the critical NLS. This solution is, hence, unstable to collapse as well, despite its linear 
stability. 

Most solutions were found to be unstable, and in fact, no bright soliton family solutions were found 
that are stable, although in the high-dimensional parameter space among the infinitely many solutions, 
we by no means discount the possibility that such solutions may exist. As examples we present the 
solution II. 2 for the values (a,6, m) = (-0.15,-1,0.9) as given in Eq. ()33p and the solution III for values 
{a,b,m,D) = (-0.1,1,0.4,-0.3). The solution with 6 < [solution II.2, Fig. ([6])] does exhibit collapse 
phenomenon like the unstable soliton, while the one with 6 > [solution III] does not. 

3 Solutions of 02 - 04 _ ^6 _ ^ 2 g^^^ hence CQNLS Model 

Let us consider solutions of field Eq. (fTOj) . As explained above, once these solutions are obtained, then the 
solution of the CQNLS equation are immediately obtained from here by using Eqs. ([T]), ([5]), ([6]), ^ and 
([5]). We list below five distinct solutions to the field Eq. (|10|) out of which four are periodic. 
Solution I 
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It is easily shown that 



Asn{Gx + xo,m) + B ^^^^ 



y Dsn{Gx + xq, m) + F 

is an exact solution to the field Eq. (jlOp provided the following five coupled equations are satisfied: 

mG^ {AF - BD){3AF + BD) =g^D'^- AA^ {aD'^ + bAD + cA^) , (48) 

G'^{AF - BD)[2mBF + (1 + m)AD] = 2g'^D^F - AaAD^AF + BD) 

-2bA^{AF + 3BD) -8cA^B , (49) 

(1 + m)G'^{AF - BDf = -Qg'^D'^F'^ + Aa{A^F^ + B^D"^ + AABDF) 

+l2hAB{AF + BD) + 2AcA^B^ , (50) 

G'^iAF - BD)[2AD + (1 + m)BF] = -2g^DF^ + AaBF{AF + BD) 

+2bB'^{3AF + BD) + 8cAB^ , (51) 

- G^{AF - BD){AF + 3BD) = g^F^ - AB^iaF^ + bBF + cB"^) . (52) 

Thus we have five coupled equations relating the five parameters A, B , D,G, g (note that we can always 
remove one of the parameter from the ansatz ()47p ). In the special case of m = 1, the solution ()47p goes 
over to the hyperbolic soliton solution 



Atanh(Ga; + xq) + B 



y D tanh(Gx + xq) + F 
while in the limit m = it goes over to the trigonometric solution 



(53) 



ylsin(Gx + xo) + B 



y Dsm{Gx + Xq) + F ' 
In the special case of -D = 0, -P = 1, we obtain a simpler solution 



(54) 



(j) = -^AsniGx + xo,m) + B , (55) 

provided 



^2 962-32|a||c| ,o SmG"^ „ 3b 



a<0, c<0, G'= ' " ' , A' = ——, B = —, (56) 

8|c|(l + m) 4|c| 8|c| 
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while g for this solution is given by 



^2 = 4^2^^ + hB + cB^) - A^G'^ 




-{U\a\\c\fm - 96^(26m - 5 - 5m^)] . 



(57) 



It is easily checked that so long as g ^ such a solution does not exist for m = 1, even though it exists 
for a range of values of m < 1. Also note that the solution (|55p does not exist in the trigonometric limit 
of m = unless c = 0. 
Solution II 

Yet another solution to the field Eq. (|10p is 




(58) 



provided the following five coupled equations are satisfied: 



mG'^{AF - BD){3AF + BD) 



g'^D^ + AA^{aD'^ + hAD + cA^) 



(59) 



G'^{AF - BD)[2mBF + (2m - l)^^] = -Ig'^D^F 



+4:aAD{AF + BD) + 2bA^{AF + 3BD) + 8cA^B 



(60) 



(2m - 1)G'^{AF - BDf = 6g^D^F^ - AaiA^F'^ + B'^D^ + 4ABDF) 



12bAB{AF + BD) - 2AcA^B^ 



(61) 



G'^iAF - BD)[2AD{1 - m) - (2m - 1)BF] = -2g^DF^ 



+4aBD{AF + BD) + 2bB^{3AF + BD) + 8cAB^ , 



(62) 



(1 - m)G^{AF - BD){AF + 3BD) = g^F'^ - 4:B^{aF'^ + bBF + cB"^) . 



(63) 



In the special case of m = 1, the solution ()58p goes over to the hyperbolic soliton solution 




(64) 
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while in the hmit m = it goes over to the trigonometric solution 



Acos{Gx + xo) + B 



y Dcos{Gx + xo) + F ' 
In the special case of -D = 0, F = 1, we obtain a simpler solution 



(65) 



(f) = y/Acn{Gx + xo,m) + B , (66) 

provided 

8c(2m - 1) ' 4c 8c ' ^ ^ 

while g for this solution is given by 

= AB^ia + bB + cB^) - (1 - m)A^G^ 

-{QAacfmil - m) - 96^(5 + 16m - 16m^)] . (68) 

Thus unlike solution ()47p . such a solution always exists in the hyperbolic limit so long as 32ac < 9b^ < 
(192/5)ac. Further, it also exists for a range of values of m, though it does not exist in the m = limit, 
unless c = 0. 

In the special case of m = 1/2, we obtain a one-parameter family of solutions of the form (j66p . In 
particular, in that case it follows from Eq. (j67|) that 

962 = 32ac, = — , B = -- , g' = — . 69 

8c 8c 3c 3 

It is worth noting that while A, G are arbitrary with their ratio being fixed, the only constraint that we 

have is that vl^ < ^ so that g is nonzero and real. 

Unlike the sn case, in this case solutions exist even in case A = (and assuming B = 1 without any 

loss of generality) . In particular, 

(p = I , (70) 

y/ Dcn{Gx + xo,m) + F 

is an exact solution to Eq. ([TOD provided 

mG^ = g^D^ , 3(1 - m)D^G^ = g'^F^ - A{aF^ + bF + c), 

(2m - 1)G^ = 2(3^2^2 - 2a) , Ag^F^ = b + AaF . (71) 
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Note that such a solution exists in the hyperbohc hmit (m = 1) but not in the trigonometric hmit {r. 
Solution III 

Yet another solution to the field Eq. (fTUj) is 



' Adn{Gx + xo,m) + B 



y Ddn{Gx + Xq, m) + F ^ 
provided the following five coupled equations are satisfied: 

- G^{AF - BD){3AF + BD) = g'^D^ - AA^iaD^ + bAD + cA^) , 

G'^iAF - BD)[2BF + (2 - m)AD] = -2g'^D^F + 4.aAD{AF + BD) 
+2bA'^{AF + WD) + ScA^B , 

(2 - m)G^{AF - BD)"^ = Gg^D^F"^ - AaiA^F"^ + B'^D^ + AABDF) 
-12bAB{AF + BD) - 2AcA^B'^ , 

G'^{AF - BD)[2AD{1 - m) - (2 - m)BF] = -2g'^DF^ 
+AaBD{AD + BF) + 2bB'^{3AF + BD) + ScAB^ , 

(1 - m)G^{AF - BD){AF + 3BD) = g^F^ - AB^{aF^ + bBF + cB"^) . 

In the special case of m = 1, the solution ()72p goes over to the hyperbolic soliton solution ([6; 
In the special case of D = 0, F = 1, we obtain a simpler solution 

(p = a/ Adn{Gx + Xq, m) + B , 

provided 

. 9b^-32ac . 3G^ 3b 
8c(2 — m) 4c 8c 

while g for this solution is given by 

= 4^2(0 + bB + cB^) + (1 - m)^^^^ 
3 

■[192b^ac{m'^ + 8m - 8) 



[32(2 - m)]2c3 ' 
+(64ac)2(l -m) + 96*^(16 - 16m - Sm^)] . 
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Like the cn case, in this case too solutions exist in case ^4 = (and assuming B = 1 without any loss 
of generality). In particular, 

^ (81) 



y/Ddn{Gx + xo, m) + F ' 
is an exact solution to Eq. (llOh provided 

= g'^D^ , 3(1 - m)D'^G^ = -g'^F^ + 4{aF^ + bF + c) , 
(2 - m)G^ = 2(35^^2 - 2a) , 4g^F^ = b + 4aF . (82) 

Note that this solution is valid for all values of m including the hyperbolic limit of m = 1. 

However, unlike the cn case, in this case solution exists even in case B = (and assuming A = 1 
without any loss of generality) . In particular 



' dn(Gx + Xo, m) 
Ddn{Gx + Xo, m) + F ' 



(83) 



is an exact solution to Eq. pO|) provided 



(1 - m)G^ = g^F^ , 3F^G^ = -g'^D^ + 4:{aD^ + bD + c) , 

{2 - m)G^ = 2{3g^ D"^ - 2a) , Ag^D^ = b + 4aD. (84) 

Note that such a solution does not exist in the hyperbolic limit so long as g ^ 0. 
Solution IV 
It is easily shown that 



(85) 



Asn'^{Gx + xojm) + B 
^ ~ ^ Dsn^{Gx + xo,m) + F' 

is an exact solution to the field Eq. provided the following five coupled equations are satisfied: 

4mG'^AD{AF - BD) = -g^D^ + AA'^{aD'^ + bAD + cA^) , (86) 

G'^iAF - BD)[2mAF + 2(1 + m)AD] = g'^D^F - 2aAD{AF + BD) 

-bA^iAF + 3BD) - AcA^B , (87) 
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2G^{AF - BD)[3mBF - (1 + m)AF + (1 + m)BD - 3AD] 
= 3g^D^F^ - 2a{A^F'^ + B^D^ + AABDF) 
-QhAB{AF + BD) - UcA^B"^ , 

2G'^{AF - BD)B[D + (1 + m)F] = -g'^DF^ + 2aBF{AF + BD) 
+bB^{3AF + BD) + AcAB^ , 

4G'^{AF - BD)BF = g^F^ - AB'^{aF^ + bBF + cB'^) . 



Thus we have five coupled equations relating the five parameters A, B, D, G, g. In the special 
the solution ([85]) goes over to the hyperbolic soliton solution 



while in the limit m = it goes over to the trigonometric solution 



Solution V 

The field Eq. (jlOp also has a remarkable nonperiodic solution given by 



g^D^ = 4A^{aD^ + bAD + cA^) , 
g'^D^F = 2aAD{AF + BD) + bA^iAF + 3BD) + AcA^B , 

6AD{AF - BD) + Sg'^D^F'^ = 2a{A'^F^ + B'^D^ + AABDF) 
+6bAB{AF + BD) + UcA^B^ , 

2BD{AF - BD) + g'^DF^ = 2aBF{AF + BD) 
+bB^{3AF + BD) + AcAB^ , 

c,2i?4 ^ 4BD{AF - BD) + AB'^{aF^ + bBF + cB"^) . 






provided the following five relations are satisfied 
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3.1 Stability of nonlinear-phase- modulated solutions 

We now briefly discuss the stability of the nonhnear-phase-modulated solutions presented in this section. 
We emphasize that dynamical instability is unaffected by a non-trivial linear phase factor, corresponding 
to a traveling frame x ^ x — vt as described in Sec. [TJ However, it is not clear a priori what the physical 
relevance may be of a solution with nonlinear phase, such as thosepresented here. Also, none of the 
solutions presented here are related in any way to the stable family of solutions from the preceding section. 
In fact, only Solution I could be considered among the dark soliton family, but it does not exist for m = 1, 
and is not similar in form to the stable solution from the previous section. Additionally, the question of 
finding these solutions numerically is significantly more challenging than the trivial phase solutions, at 
least for the periodic amplitude ones, owing to the necessity of finding suitable parameters which satisfy 
systems of coupled nonlinear equations and lead to admissible (real-valued) and nontrivial solutions. We 
do not expect these solutions to be more stable than their trivial phase counterparts from the previous 
section. Therefore, while the solutions presented herein are valid and interesting in their own right, we do 
not investigate their stability in detail here. Nonetheless, for illustrative purposes we briefly explain the 
methodology and then present two solutions and their stability, one at the hyperbolic limit, m = 1, and 
one elliptic solution. 

The criteria in this section for the 10 parameters involved in the equation and corresponding solution 
consist, in general, of 5 coupled nonlinear equations, for which one can choose 5 parameters and solve for 
the remaining 5 with Newton's method, for instance. In this way, one often flnds that A/B=D/F, which 
is a trivial solution, or that the parameters either are imaginary or lead to imaginary solutions. 

Once the solution cp to Eq. [TU]is found, the solution u = (pexp {ig/2ri} must be found in a numerical 
domain in order to examine its stability. This leads to the additional resonance condition that there exist 
integers n and k such that 

rnr/g = kP{m,G) (99) 

where P{m,G) = 4:K{m)/G is the period of the elliptic function (j) and K{m) is the complete elliptic 
integral of the flrst kind. The domain has to be then truncated at the period of r/, and of (j) iii case 
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Figure 7: (Color Online) The top panels depict solution and phase and the bottom panels show the linear 
spectra. The left-hand solution is the elliptic solution given in Eq. (j66p with (a, b, c, m, g, A, B, D, F, G) = 
(0.5077, -0.9501, 0.5, 0.5, 0.4145, 0.1, 0.7126, 0, 1, 0.1155), while the right-hand solution is for the hyperbolic 
limit with {a,b,c,m,g,A,B,D,F,G) = (1,-2,1,1,0.375,0.6124,0.75,0,1,0.7071). 

m ^ 1. If m = 1, we choose one period of r] as the domain, i.e. {L;r] = ir/g}. For m 7^ 1, we determine 
(n*,k*) = argmin{(n. A:); nvr^ = kP{m,G)} and then truncate the domain at {L;r] = mr/g}. 

As displayed in Fig. [71 the stability does not differ significantly qualitatively from the solutions with 
trivial phase (compare with the solutions in Figs. l4|6l top). except the hyperbolic solution is now periodic 
in phase, and, hence, has a continuum of unstable eigendirections. The dynamics of both these solutions, 
when perturbed with additive noise, results in collapse, as in the case of trivial phase. 
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4 Conclusions 



To conclude, a large array of exact analytical solutions to various field theories were found leading directly 
to various exact solutions of cubic-quintic NLS and their stability and dynamics were studied numerically. 
The majority of these are oscillatory solutions which are unstable, however, some are found to be linearly 
stable, such as some elliptic solutions from the dark-soliton family solution II for appropriate parameter 
regimes. Also many hyperbolic dark-solitons were found to be stable and quite robust over a range of 
parameter values. The hyperbolic bright soliton solutions were found to be unstable to collapse when the 
cubic term is defocusing or nonexistent and the quintic term is focusing, while they were found to be always 
stable when the cubic term is focusing, for both focusing and defocusing quintic term (at least within the 
parametric regions examined herein). Some periodic solutions with defocusing and focusing cubic and 
quintic nonlinear terms, respectively, also displayed collapse phenomenon, while ones with a focusing cubic 
term did not. 

A Solutions of 0^-0^-0^ and hence Quadratic- Cubic NLS Model 

Let us consider solutions of field Eq. (fT2|) . As explained above, once these solutions are obtained, then the 
solution of the quadratic-cubic equation are immediately obtained from here by using Eqs. ([T]), ([2]) and 
(jlip . We list below eight distinct periodic solutions to the field Eq. ()12p . In each case, we also mention 
the values of the parameters a, b, c, in particular, if they are positive or negative. However, instead of Eq. 
()12p we use slightly different form given by 



This is done so that one can easily pick up the solutions recently obtained by us in a related paper on 
coupled (jP' -(jr" -(j)'^ field theory [9]. 
Solution I 

It is easily shown [9| that 



(t)xx + 2a(/> + 36(/>^ + 4c(/)^ = . 



(100) 




(101) 
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is an exact solution to the field Eq. (jlOOp provided 

{l + m)B'^ = -a, 2{1 + m)cA^ = ma , b = -4cF , l? = Aac. (102) 

Thus this solution is valid provided a < 0, c < while b could be positive or negative. 
In the limit m = 1 the periodic solution (jlOip goes over to the dark soliton solution 

(t) = F + A tanh(53; + xq) , (103) 

provided 

25^ = -a , AcA^ = a , b = -AcF , b^ = Aac . (104) 

Thus the dark soliton solution exists to field Eq. (jlOOp provided a < 0, c < while b could be positive or 
negative. 

Solution II 

It is easily shown that 

(/> = F + ^cn(Sx + xo,m), (105) 
is an exact solution to the field Eq. (jlOOp provided 

(2m - \)B^ = a , 2cA'^ = mB^ , b^ = Aac,b = -4cF . (106) 

Thus this solution is valid provided a > 0, c > while b could be positive or negative. Note also that this 
solution exists only if m > 1/2. 

In the limit m = 1 the periodic solution (jlOSp goes over to the bright soliton solution 

4> = F + Asech{Bx + xq) , (107) 

provided 

B"^ = a, 2cA^ = B^ b^ = Aac, b = -AcF . (108) 

Thus the bright soliton solution exists to field Eq. (jlOOp provided a > 0, c > while b could be either 
positive or negative. 
Solution III 
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It is easily shown f9] that 

(t> = F + Adn{Bx + xo,m), (109) 
is an exact solution to the field Eq. (jlOOp provided 

(2m- 1)5^ = 0, 2c^2 = b2, b'^ = 4:ac, b = -AcF . (110) 

Thus this solution is valid provided a > 0, c > while b could be either positive or negative. Note also 
that unlike solution II, this solution is valid for any m(0 <m< 1). 

In the limit m = 1 the periodic solution (|109p goes over to the bright soliton solution (|107p and hence 
satisfy constraints given by Eq. (jlOSp . 

Solution IV 

It is easily shown [9j that 

c^ = F+ -4sn(Z).. + xo,m) ^^^^^ 
1 + Bdn[Dx + xo, m) 

is an exact solution to the field Eq. (llOOp provided 

B = l, (2 - m)L»2 = -2a, ^{2 - m)cA'^ = m'^a , b^ = 4.ac, 5 = -4cF . (112) 

Thus this solution is valid provided a < 0, c < while h could be either positive or negative. 
In the limit m = 1 the periodic solution (jllip goes over to the dark soliton solution 

<t^ = F+ -4tanh(^x + xo) 
^ 1 + Bsech{Dx + xq) 

provided 

B = l, D^ = -2a, AcA^ = a, b^ = 4ac, b = -4cF . (114) 

Thus the dark soliton solution exists to field Eq. (jlUOp provided a < 0, c < while b could have either sign. 
Solution V 

It is easily shown [9] that 

^cn(L>x + xo,m) 
V 1 — m + Ban[L>x + xq, m) 

is an exact solution to the field Eq. (jlOOp provided the parameters satisfy the same constraints as given 

by Eq. (jll2p . Thus this solution is valid provided a < 0, c < while b could have either sign. 
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In the limit m = 1 the periodic solution (jllSp goes over to a constant, i.e. 

(t> = F + A, (116) 

where F, A are as given by Eq. (jll4p . 
Solution IV 

It is easily shown |9j that 

Atanh{Dx + xq) 

^ ^ l + Btanh(Z)x + xo) ' ^ ' 

is an exact solution to the field Eq. (jlOOp provided 

2D^ = -a, A = F(1-B), b'^ = 4ac, Fb = -{l + B)a. (118) 

Thus this solution is valid provided a < 0, c < while h could be either positive or negative. 
Solution VII 
It is easily shown that 



(119) 



^ Asech.{Dx + xq) 
1 + Bsech.{Dx + xq) ' 

is an exact solution to the field Eq. (jlOOp provided 

A = ^==, D^ = -2a, B= . (120) 

Thus this solution is valid provided a < while b, c could have either sign. In particular, if c > then b 
can have either sign while if c < then 6 > 0. This solution can also be rewritten as 

-2a 



b + Vb"^ — 4accosh[v2a(x + xq)] 

Solution VIII 

It is easily shown [9j that 



(121) 



^ Asedi{Dx + xo) 

1 + Bsech{Dx + xq) ' ^ ' 



is also an exact solution to the field Eq. (jlOOp provided 



^ = -J^±^, D^ = 3bF + Aa, B = -4±^. (123) 
^/b^ + 2bcF VWT2bcF 
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Thus this solution is vahd provided 6 > 0, c < while a > (< or =)0 depending on if m < (> or =)l/5. 
Solution IX 
It is easily shown that 

Ax'^ + B , , 

is an exact solution to the field Eq. pOOp provided 

AE = -3BD, ;g=-y, 9b^ = ^2ac, ■^ = -^- (125) 
Thus this solution is only valid if a > 0, c > while b could be either positive or negative. 



B Solutions of 0^-0^ and hence Standard (Cubic) NLS Model 

Finally, for completeness, we consider solutions of field Eq. (I12p . in case 6 = 0. As explained above, 
once these solutions are obtained, then the solution of the standard (cubic) NLS equation are immediately 
obtained from here by using Eqs. ([T]), ([2]) and (jllh . We list below three distinct periodic solutions to the 
field Eq. (|12p with 6 = 0. In each case, we also mention the values of the parameters a and c, in particular, 
if they are positive or negative. However, instead of Eq. (I12p (with 6 = 0) we use slightly different form 
given by 

(l)a;x + 2a(l) + Accj)^ = . (126) 

This is done so that one can easily pick up the solutions recently obtained by us in a related paper on 
coupled (/)^ — (j)^ field theory |lOj . 
Solution I 

It is easily shown [TO] that 

4> = Asn{Bx + xo,m), (127) 
is an exact solution to the field Eq. (jl26p provided 

{l + m)B^ = 2a, {1 + m)cA'^ = -ma . (128) 

Thus this solution is valid provided a > 0, c < 0. 
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In the limit m = 1 the periodic solution (jl27p goes over to the dark soliton solution 

(j) = Atanh{Bx + xo) , (129) 

provided 

B^ = a, 2cA^ = -a. (130) 

Thus the dark soliton solution exists to field Eq. ()126p provided a > 0, c < 0. 
Solution II 

It is easily shown |10j that 

(f) = AcTi{Bx + XQ,m) , (131) 
is an exact solution to the field Eq. (jl26p provided 

(2m- 1)^2 = -2a, {2m - l)cA^ = -ma . (132) 

Thus this solution is valid provided c > while a < (> or =)0 depending on if m > (< or =)l/2. 
In the limit m = 1 the periodic solution (|13ip goes over to the bright soliton solution 

(j) = Asech{Bx + Xo) , (133) 

provided 

5^ = -2a , cA^ = -a . (134) 

Thus the bright soliton solution exists to field Eq. ()126p provided a < 0, c > 0. 
Solution III 
It is easily shown |10j that 

(f) = Adii{Bx + xo,m) , (135) 
is an exact solution to the field Eq. (jl26p provided 

(2-m)B2 = -2a, {2-m)cA^ = -a. (136) 

Thus this solution is valid provided a < 0, c > 0. 

In the limit m = 1 the periodic solution (|135p goes over to the bright soliton solution (|133p and hence 
satisfy constraints given by Eq. (jl34p . 
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